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Abstract 

The free energy between a static quark and an antiquark is studied by using the color-singlet 
Polyakov-line correlation at finite temperature in lattice QCD with 2+1 flavors of improved Wil- 
son quarks. From the simulations on 32 3 x 12, 10, 8, 6, 4 lattices in the high temperature phase, 
based on the fixed scale approach, we find that, the heavy-quark free energies at short distance 
converge to the heavy-quark potential evaluated from the Wilson loop at zero temperature, in 
accordance with the expected insensitivity of short distance physics to the temperature. At long 
distance, the heavy-quark free energies approach to twice the single-quark free energies, imply- 
ing that the interaction between heavy quarks is screened. The Debye screening mass obtained 
from the long range behavior of the free energy is compared with the results of thermal pertur- 
bation theory. 



1. Introduction and the method 

Heavy-ion collision experiments running at RHIC have provided various information about 
the new state of matter, the quark-gluon plasma (QGP). To remove theoretical uncertainties in 
understanding the nature of QGP, first principle calculations by lattice QCD are indispensable. 
So far, most lattice studies at finite temperature have been performed using staggered-type quark 
actions with the fourth-root trick, whose theoretical basis is not fully established yet. Thus, a 
crosscheck with other actions such as the Wilson-type quark actions are important to control and 
estimate systematic errors due to lattice discretization. 

We investigate thermodynamic properties of QGP using the improved Wilson quark action. 
Adopting the fixed scale approach developed in OJJ], we perform finite-temperature simulations 
J2], while the corresponding zero-temperature configurations are taken from the results of the 
CP-PACS and JLQCD Collaborations fll. 

At T = 0, interaction between a static quark and an antiquark can be studied by the heavy- 
quark potential evaluated from the Wilson loop operator. The resulting potential takes the 
Coulomb form at short distances due to perturbative gluon exchange, while it takes the linear 
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form at long distances due to confinement: 



V(r) = -—+(Tr+Vo. (1) 
r 

For T > 0, inter-quark interaction may be studied by the heavy-quark free energy F l (j, T) eval- 
uated from a Polyakov-line correlation in the color-singlet channel with Coulomb gauge fixing 
ft: 

F\r, T) = -7-ln(TrO+(x)0(y)>, r = |x- y|, (2) 

where Q(x) = YljL\ U4(x,t) with t/4(x, t) being the link variable in the temporal direction. At 
zero temperature, we expect F (r, T = 0) = V(r), while at high temperature, we can extract the 
Debye screening mass ma from 

F\ r , T) = -^D. e -' nDiT)r + 2F Q . (3) 
r 

In the conventional fixed N, approach where T is varied by changing the lattice spacing a, 
V(r) and F\r, T) receive different renormalization at each T, and thus are usually adjusted by 
hand such that V(r) and F 1 (r, T) coincide with each other at a short distance, assuming that 
the short distance properties are insensitive to the temperature. On the other hand, in the fixed 
scale approach, the temperature T - (A^a) _I is varied by changing the temporal lattice size N, 
at fixed a QJL£] . Because the spatial volume and the renormalization factors are common to all 
temperatures, we can directly compare the free-energies at different temperatures without any 
adjustment. We show below that F \r, T)'s for different T approach to V(r) at short distances, 
which proves the expected insensitivity of the short distance physics to the temperature. 



2. Results of the lattice simulations 

We employ the renormalization group improved gluon action and 2 + 1 flavors of nonpertur- 
batively (9(a)-improved Wilson quark actions. Zero-temperature configurations are given by the 
CP-PACS and JLQCD Collaborations H at tn K /m p = 0.6337(38) and m K /m K . = 0.7377(28). 
Finite temperature simulations with the same parameters are performed on 32 3 x N, lattices with 
N, = 12, 10, 8, 6 and 4, which correspond to T ~ 200-700 MeV Q]. The absolute scale is 
estimated from the Sommer parameter, ro = 0.5 fm. 

Figure Q] shows the results of the heavy-quark potential V(r) at T — and the heavy-quark 
free energies F l {r, T) at various temperatures as functions of r. At T = 0, V(r) shows Coulomb- 
like and linear-like behaviors at short and long distances, respectively. A fit with Eq. (Q3, as 
shown by the dashed line in Fig.Q] gives oq = 0.441 and yjtr = 0.434 GeV. 

For T > 0, we note that the heavy-quark free energies F 1 (r, T) at all temperatures converge to 
V(r) at short distances. This means that the short distance physics is insensitive to temperature. 
As stressed above, unlike the case of the conventional fixed-A^ ; approach in which this insensitiv- 
ity is assumed and used to adjust the constant terms of F l {r, T), our fixed scale approach enabled 
us to directly confirm this theoretical expectation. 

At large r, F l {r, T) departs from V(r) and eventually becomes flat due to Debye screening. In 
Fig.Q] the asymptotic values of F l (r, T) at long distance are also compared with 2Fq denoted by 
the arrows: Here the thermal average of a single Polyakov-line is defined as Fq = -rin(TrO). 
We find that F l (r, T) converges to 2Fq quite accurately at long distances. 
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Figure 1 : Heavy-quark free energies at various temperatures. The heavy-quark potential at T = was calculated by the 
CP-PACS and JLQCD Collaborations d. 



In order to extract the screening mass mo, we fit F x (r, T) by the screened Coulomb form, 
Eq. (01. The fit range is chosen to be 0.38 fm < r < 0.57 fm. The left panel of Fig. [2] shows 
mo(T)/T which does not have strong dependence on T. To make a quantitative comparison to 
the result of the thermal perturbation theory, we define the 2-loop running coupling by g^in) = 
y6 ln(ju/A^) 2 + & In [ln^/A^) 2 ] with the QCD scale parameter A^_ =3 = 260 MeV fl, where 
we assume a degenerated Nf = 3 case. Then, the Debye mass in the leading-order (LO) thermal 
perturbation theory is given by m^°(7')/7' = Jl + 7V//6 g2\(T) neglecting the quark mass effects. 
A formula in the next-to-leading-order (NLO) is also available from the resummed hard thermal 
loop calculation J6J]: 



1 + -f 8n(T) 
o 



i+S2iCO^-' ' 



2n V 1 +7V//6 



2m h0 1 ~\ 



+ o(g 2 ) 



(4) 



where m mag (T) = C m g 2 (T)T is the magnetic screening mass. Since the factor C m cannot be 
determined in perturbation theory due to the infrared problem, we adopt C m - 0.482 calculated 
in a quenched lattice simulation [7] as a typical value. In Fig. [2jleft), the dashed lines represent 
the LO results for Nf = 3 and the bold lines represent the NLO results for Nf = 3, for a range of 
the renormalization point /v = nT-'inT . We find that the LO Debye mass does not reproduce the 
lattice data in magnitude, while the NLO Debye mass is approximately 50 % larger than the LO 
Debye mass and shows a better agreement with the lattice data. 

Finally we study the flavor dependence of mo- Fig. |2jright) shows mo in Nf = 2 + 1 QCD 
(this work), that in Nf — 2 QCD with an improved Wilson-quark action at m K /m p = 0.65 JUl, 
and that in the quenched QCD (Nf — 0) QJJ] . Arrows on the horizontal axis indicate critical 
temperatures for Nf - 2 (T c ~ 186 MeV at m n /m p = 0.65) and Nf - (T c ~ 290 MeV). We 
find that mo for Nf = 2 + 1 is comparable to that for Nf = 2, whereas it is larger than that for 
Nf - 0. A similar result was obtained with a staggered-type quark action (9J]. 

3. Summary 

We studied the free energy between a static quark and an antiquark at finite temperature in 
lattice QCD with 2+1 flavors of improved Wilson quarks on 32 3 x 12-4 lattices in the high 
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Figure 2: Left: Results of the Debye mass for Nj ■ = 2 + 1 obtained from the heavy-quark free energy with those of the 
leading order (LO) and next-to- leading order (NLO) thermal perturbation theory for Nf = 3. Right: Comparison among 
the Debye masses for iVy = 2+ l,iVy = 2(^| and Nf = extracted from the heavy-quark free energy on the lattice. 
Arrows in the right panel indicate critical temperatures in Nf = 2 and cases. 



temperature phase. We adopted the fixed scale approach which enables us to compare the free- 
energies at different temperatures directly without any adjustment of the overall constant. At 
short distances, the heavy-quark free energies, evaluated from the Polyakov-line correlations in 
the color-singlet channel, show universal Coulomb-like behavior common to that of the heavy- 
quark potential at zero temperature. This is in accordance with the expected insensitivity of short 
distance physics to the temperature. At long distances, the heavy-quark free energies approach 
to twice the single-quark free energies calculated from the thermal average of a Polyakov-line. 
Also, we extracted the Debye screening mass mo(T) from the heavy-quark free energy and found 
that the next-to-leading order thermal perturbation theory is required to explain the magnitude of 
mo(T) on the lattice. Comparison to the previous results at Nf = 2 and Nf = 0, shows that the 
dynamical light quarks have sizable effects on the value of mo(T). 
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